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Abstract 

The solution to fine tuning is one of the principal motivations for Beyond the 
. Standard Model (BSM) Studies. However constraints on new physics indicate that 

I many of these BSM models are also fine tuned (although to a much lesser extent). 

' To compare these BSM models it is essential that we have a reliable, quantitative 

measure of tuning. We review the measures of tuning used in the literature and 
Q , propose an alternative measure. We apply this measure to several toy models and 

I the Minimal Supersymmetric Standard Model. 

OO 

I Introduction 

Fine tuning appears in many areas of particle physics and cosmology, such as the Standard 
5J^| Model (SM) Hierarchy Problem and the Cosmological Constant Problem. These problems 
^ • imply that the the universe we live in is a very atypical scenario of the theories we use to 
describe it. The contortion required to reproduce observation makes such theories seem 
^ . unnatural, motivating many studies of Beyond the Standard Model (BSM) physics. 

However many of the models constructed to solve fine tuning, also exhibit some degree 
of tuning themselves. In the absence of data, while we await the LHC, naturalness is 
CN I used to compare models and judge their viability. Great importance has been attached 
IT) • to small differences in the levels of tuning when comparing models, so it is important 

that naturalness and fine tuning are rigorously understood and measured accurately. 
Q ! For example the Hierarchy Problem is one of the fundamental motivations of low 

energy supersymmetry (SUSY) (for a review see Ref. [1]). If the SM is an effective theory, 
valid up to the Planck scale, then the inclusion of supersymmetric partners for every SM 



Oh 



X 



■ particle leads to the cancellation of quadratic divergences in the loop corrections to the 



Higgs mass. This removes the need for fine tuning of O{10^^) between the tree-level mass 
parameter and the Planck Mass, allowing the Higgs boson to be naturally light. 

Unfortunately current limits on superpartner masses may imply fine tuning in the 
most studied model, the Minimal Supersymmetric Standard Model (MSSM). The min- 
imisation of the Higgs potential sets the square of the Z boson mass, M|, in terms of the 
supersymmetry breaking scales. In the MSSM the tree-level expression for this is, 

,j2 _ 2(^L-^ktan^/^) 

^^^^ " tan2/5-l ' 

where tan (3 is the ratio of vacuum expectation values, fi the bilinear Higgs superpotential 
parameter, and ttlh^ and ttlh^ are the up and down type Higgs scalar masses respectively. 
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Lower bounds on the masses of the supersymmetric particles and the Higgs translate 
to lower bounds on the parameters appearing on the right hand side of Eq. ([T]). If, 
for example, one of the parameters is 1 TeV, then to cancel this contribution and give 
Mz = 91.1876±0.0021 GeV [2], another parameter (or combination of parameters) would 
have to be tuned to the order of one part in a hundred. 

Including loop corrections to Eq. ([T]) and examining the experimental constraints, one 
finds that the largest term is from corrections involving the heaviest stop. This can be 
written as |3], 



where A is the high scale at which the soft stop masses, m^^ and m^-^ , are generated from 
the supersymmetry breaking mechanism and yt is the top Yukawa coupling. A heavy 
physical stop mass (m^ > 500 GeV) is needed to provide radiative corrections to the light 
CP even Higgs mass, m/jO of the form, 

^"^ho = -j^ sm /?log j , (3) 

which are large enough to evade the LEP constraints on it's mass (> 114 GeV). So the 
Little Hierarchy Problem is really about the tension between the masses of the Z boson, 
the heaviest stop squark and the light Higgs. 

The desire to solve this "Little Hierarchy Problem" has motivated a flood of activity 
in the construction of supersymmetric models [IHTU]. There is also increased interest in 
studying alternative solutions to the SM Hierarchy problem [TTHT3] . In addition to ensur- 
ing such models satisfy phenomenological constraints it is essential that the naturalness 
is examined using a reliable, quantitative measure of tuning. 

In Ref. |il4j Barbieri and Giudice use a measure of tuning, originally proposed in Ref. |15j. 
for an observable, O, with respect to a parameter, pi, 

A large value of AsciPi) implies that a small change in the parameter results in a large 
change in the observable, so the parameters must be carefully "tuned" to the observed 
value. Since there is one ABoiPi) per parameter, they define the largest of these values 
to be the tuning for that point, 

Abg = max({ABG(Pi)}). (5) 

They then make the aesthetic choice that a tuning, A^g > 10 is fine tuned. 

This measure has been used extensively in the literature to quantify tuning in the 
MSSM p!6ti26] and to examine tuning in other models and theories e.g. [3]. [27143T] . How- 
ever other measures have also been proposed and used in the literature. 

Motivated by global sensitivity, which will be discussed in the next section, Anderson 
and Castano [321435] propose that tuning should be measured with, 

AAciPi) = ^ (6) 

A EG [Pi) 
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where they choose the "average" sensitivity, ABciPi), not to be the mean, but instead 
defined by, 

pf{p)Jdp 

where f(p) is the probabihty distribution of parameter p. Individual A^cfe) are com- 
bined in the same manner as the individual /^BciPi) were. 



Aac = max({AAc(Pi)}). (8) 

There is some dispute within the literature as to whether or not Eq. ([5]) is the best way 
of choosing a final tuning value from the set {ABciPi)}- In [TT] . [361439] the individual 
ABciPi) are be combined as if uncorrected. 



^E= ^J2^BGiP^)■ (9) 

Several other measures have been proposed [i0lfl3] . but will not be discussed here. 

In Section [IT] we detail some limitations of the traditional measure of tuning, Abgj used 
for this Little Hierarchy Problem. We then describe our fundamental notion of tuning, 
and how this principle can be applied to construct quantitative measures of tuning in 
Section IIIII This leads us to present a new tuning measure in Section IIVI which is also a 
generalisation of the traditional measure that overcomes the limitations outlined earlier. 
This measure is applied to several toy models in Section |V] to demonstrate how it works 
and compare the results it produces with those from other measures. Finally in Section 
IVII we apply our measure to the Little Hierarchy Problem for a selection of Minimal 
SUperGRAvity (MSUGRA) inspired points. 

II Limitations of the Traditional Measure 

Despite the wide use of A bg it has several limitations which may obscure the true picture 
of tuning: 

• variations in each parameter are considered separately; 

• only one observable is considered in the tuning measure, but there may be tunings 
in several observables; 

• it does not take account of global sensitivity; 

• only infinitesimal variations in the parameters are considered; 

• there is an implicit assumption that the parameters come from uniform probability 
distributions. 

Tuning is really concerned with how the parameters are combined to produce an unnatural 
result. If one measures tunings for each parameter individually, there is no clear guide 
how to combine these tunings to quantify how unnatural this cancellation is. This has 
led to two alternative approaches in the literature, Eq. ([5]) and Eq. ([9]); the only way to 
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determine if either Abg ot Ae combines sensitivities correctly is to compare them with 
a generahsation of A^cfe) that varies all of the parameters simultaneously. 

Secondly, some theories may contain significant tunings in more than one observable. 
We want to know how can these tunings be combined to provide a single measure. For 
example it is reported in Refs. [11113], and more recently in Refs. |46|ITF] . that the MSSM 
also requires tuning in the relic density of the dark matter (p). To measure the tuning 
for some particular set, S' = {M'^,p'}, of these observables we should determine how 
atypical predictions like S' are in the theory. There are four classes of scenario which 
are significant: the first where both Mz and p are similar to their value in S'; two more 
classes where only one of Mz or p is similar to it's value in S'; one with neither observable 
similar to 5". Tunings in these two observables should be combined in a manner which 
measures how atypical scenarios in the first class are, without double counting scenarios 
which appear in the final class. Only a tuning measure which considers the observables 
simultaneously can achieve this. 

A third problem, first mentioned by Anderson and Castano [32j is that the traditional 
measure picks up global sensitivity as well as true tuning. Abg is really a measure of 
sensitivity. Consider the simple mapping / : x — x", where n ^ 1. Applying the 
traditional measure to f{x) gives Abg = Abg{x) = n. Since Abg is independent of x, we 
follow the example of p2] and term this global sensitivity. Since Abg{xi) — Abg{x2) = 
for all Xi,X2, there is no relative sensitivity between points in the parameter space. 

If we use A^G' as our tuning measure then /(x) appears fine tuned throughout the 
entire parameter space. This contrasts with our fundamental notion of tuning being a 
measure of how atypical a scenario is. A true measure of tuning should only be greater 
than one when there is relative sensitivity between different points in the parameter space. 

Another concern is that A bg only considers infinitesimal variations in the parameters. 
Since MSSM observables are complicated functions of many parameters, it is reasonable 
to expect some complicated distribution of the observables about that parameter space. 
There may be locations where some observables are stable (unstable) locally, but unstable 
(stable) over finite variations. 

Finally, there is also an implicit assumption that all values of the parameters in the 
effective softly broken Lagrangian Csusy are equally likely. However they have been 
written down in ignorance of the high-scale theory, and may not match the parameters 
in, for example, the Grand Unified Theory (GUT) Lagrangian, Cgut- Any non-trivial 
relation between these different sets of parameters may alleviate or exacerbate the fine 
tuning problem. 

While some of the alternative measures in the literature are motivated by one of these 
issues, no proposed measure fully addresses all of them. 



Ill Constructing Tuning Measures 

A physical theory is fine tuned when generic scenarios of the theory predict very dif- 
ferent physics to that which is observed. For the theory to agree with observation the 
parameters must be adjusted very carefully to lie in an extremely narrow range of values. 
Insisting that the physics described by the theory is similar to that observed, shrinks the 
acceptable volume of parameter space. When in this tiny volume even small adjustments 
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to the parameters will dramatically change the physics predicted, so fine tuning may 
also be characterised by instability. It is this instability which the traditional measure is 
exploiting. 

Instead we wish to construct a tuning measure which determines how rare or atypical 
certain physical scenarios are. The most direct way to do this is to compare the volume 
of parameter space, G, that is similar to some given scenario with the typical volume, T, 
of parameter space formed by scenarios which are similar to each other. 

If all the parameters {pi} are drawn from a uniform probability distribution then the 
probability of obtaining a scenario in G is, G/V, where V is the volume of parameter 
space formed by all possible parameter choices. Similarly T/V gives the probability of 
obtaining a scenario in volume T. We may then define tuning as A = T/G, to quantify 
the relative improbability of scenarios similar to our given scenario in comparison to the 
typical probability. 

To place this within a quantitative framework we must define what we mean by 
"similar" and "typical". This will be dealt with later. First, though, consider the toy 
example presented in Fig. [U showing an observable, O, which depends on a parameter, x. 
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Figure I: A toy example with an observable, O, which depends on a parameter x. 

Here there are four clearly distinct groups of observable scenario (O = Oi, O2, O3, O4) and 
"similar" can be replaced with equal. Given one of these groups of scenarios, = 0,, 
the volume G is the length (one dimensional volume) of parameter space with O = Oi. 
For example, for O4 we have G = X2 — Xi. Next we must define our "typical" volume, T 
formed by these distinct groups of scenarios. In this simple example an obvious choice is 
to define T as the mean volume (length) of parameter space formed by scenarios in the 
same group. So T = |(xi + {x2 — xi) + (xa — X2) + (a;4 — X3)) = ^. The tuning required 
to get O = O4 is then A = ^^^t^^-^ , which conforms to our intuitive expectation. 

In more realistic examples the definitions of "similar" and "typical" will not be so 
trivial. The definitions must be chosen to fit the type of problem one is considering. In 
the simple example given above the problem was that scenarios where = 04 occupied 
a smaller proportion of the parameter space than other values, O = Oi, O2, O3. 

In hierarchy problems the concern is that one (or more) observable is much smaller 
than another observable, despite depending on common parameters. The requirement 
that one observable is large forces the theory into a region of parameter space where 
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generic points also predict a large value for the second observable(s). 

So "similar" must be related to the size of the observables. For example, one might 
consider "similar" to observable to mean observables "of the same order" as 0[. A 
sensible definition of G is then the volume of parameter space where < §7 < 10, for 
all observables Oj. However it is not clear that this is more appropriate than some other 
choice such as | < ^ < 2. So generally G can be defined by a class of parameter space 

i 

volumes formed from dimensionless variations in the observables a < % < h. Different 

i 

values of a and h quantify different definitions of "similar" and are therefore different 
fine-tuning questions. In comparison, the one dimensional measure /S.bg is a ratio of 
infinitesimal lengths, so implicitly adopts the choice a, b —>■ 1. One can imagine cases 
where this would be a bad choice (e.g. an observable which oscillates quickly when the 
parameter is varied), so care must be taken to choose a and b sensibly (i.e. ask the correct 
question) . 

When a large hierarchy between observables requires a large cancellation between 
parameters, as in the traditional hierarchy problem, the region of parameter space which 
can provide the correct observables (the volume G) is much smaller than one would expect 
(i.e. it is "fine-tuned"). We must compare this volume with the "typical volume" of 
parameter space, T, that one would expect if no fine-tuning were present. The remaining 
question is then, how do we define this "typical volume"? 

One might suggest that this typical volume should be the average of volumes G 
throughout the whole parameter space, (G). However, the measure would then depend 
only on how far parameters are from some hypothesised upper limits on their values. For 
example, an observable O which depends on a parameter p according to O — ap will 
display fine-tuning for small values of p if one chooses the maximum possible value of p 
to be large, even though there is no cancellation present. This is not the 'fine-tuning' 
we are trying to probe; we want to gain insight into the unnatural cancellation between 
parameters, so T must be anchored to the specific parameter point to be tested. 

We can do this by adopting the same notion of "similar" that we used to define G. 
We introduce a volume F which is formed from dimensionless variations [a, b] in the 
parameters. A comparison of F/ G at different points in the parameter space, provides a 
test of whether G"s variation is due to a simple scaling with the parameters (as described 
above for O = ap), or due to some "unnatural" effect such as fine-tuning. Consequently 
one should compare F/G with its average value over the entire space, (F/G). Reverting 
to our previous terminology, the "typical" volume which one would have expected to form 
from dimensionless variations in the parameters about {p^}, is 



IV A New Measure 

Following the above discussion and motivated by the limitations of the traditional mea- 
sure, we propose a new measure of tuning. 

We define two volumes in parameter space for every point P'{p[}. Let F be the 
volume of dimensionless variations in the parameters over some arbitrary range [a, b] , 
about point P', i.e. the volume formed by imposing a <^ <b. Similarly let G be the 
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(10) 
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volume in which dimensionless variations of the observables fall into the same range [a, b], 
i.e. the volume constrained by a < < b. Volumes F and G are illustrated for a 

two dimensional example in Fig. [Ill 



PI = fcp'i 



F 





Figure II: Left: In two dimensions the bounds placed on the parameters, a < ^ < b, 

Pi 

appear as four lines in parameter space giving the dark grey area (2d volume), F. Middle: 
Bounds on the two observables, a < < b introduce four more lines giving the 

volume G. Right: Two dimensional volumes (areas) F (dark grey) and G (light grey). 

We define an unnormalised measure of tuning with. 



A 



G' 



This is sufficient for comparing different regions of parameter space within a given model 
as the normalisation factor will be common. To compare tuning in different models we 
need to include normalisation, 

A = i^, (12) 



with, 



AG" 

F\ ^ fdpl...dpn^{{p^},{0,}) 

G 



^~\G/- fdpl..dp^ ' • ^^^^ 

Notice that this measure does not depend on experimental constraints. In naturalness 
problems such constraints should only rule out the point, P', around which we make 
variations to test fine tuning. If P' is not experimentally excluded, we should not impose 
experimental constraints on nearby points {Pi} used to probe fine tuning. Fine tuning 
quantifies how unnatural a region of parameter space is and this is a feature of the theory, 
not our experimental knowledge. 

A quantifies the restriction on parameter space. This is more in touch with our in- 
tuitive notion of tuning than the stability of the observable. Notice that with only one 
or two parameters and no global sensitivity, A bg also describes restriction of parameter 
space and yields the same results as our new measure. However it is important to recog- 
nise that A^c's ability to do this leads to its utility as a tuning measure there. Equally 
its failure to do so in many dimensions demonstrates its limitation. 



Consider fine tuning for a single observable which depends on more than one parameter. 
Even though the true tuning for any physical scenario should be described using all 
available observables, it is often useful to define individual tunings for each observable 
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separately. However, in this case, the volume G is unbounded, since a single observable 
can only constrain one combination of parameters. 

To resolve this difficulty one must either reduce the number of parameters to one or 
introduce some other bounds on G. The former reintroduces the problem of combining 
tunings for individual parameters and a better procedure is to restrict G to be within 
F. Here we are trying to pick up how much of the restriction in parameter space is due 
to this particular observable. The assumption is made that if all other observables were 
natural then they would restrict G no more than F does. Therefore we define Gqj to be 

the volume restricted by a < < b and a < < h. Tuning is then defined by, 



\Co, 

This definition is applied to obtain individual tunings in the MSSM in Section I VII 

Like Abg and /\aci A depends upon the choice of parameterisation. Since tuning is 
about the restriction of the parameter space this seems unavoidable. To examine differ- 
ent choices of parametrisation one must redefine volumes F and G in terms of the new 
parameters and normalise the tuning by taking the average over the new parameter space. 

Since much of the motivation behind developing this measure was to generalise A^^ so 
that many parameters and many observables are considered simultaneously, it is inter- 
esting to look at how the two measures are related. 

Consider a theory with one observable, y which has a linear dependence on a single 
parameter, x, with the value of that parameter being drawn from a uniform probability 
distribution. At the parameter point (xo,i/o), notice that, A^g = l^o/Uody / dAi "while we 
can see F = {h — a)Xo and G = ||(& — o)l/o, so, 

_ hxo- axo dy _ ^^^^ 
G byo - ayo dx 

Similarly Anderson and Castano's measure may be written as, 

_ Xody jdx'yjx')^^ _ J dx'yjx') ^^^^ 
yo dx Xo J dx' yo J dx' 

Now notice that (G) = gXM^z^, go, 

A.G = (17) 

In Section UTTl we pointed out the difficulty in using {G)/G as a tuning measure and 
this will be further illustrated in Section |V] when we look at results for our measure and 
Aac for a toy version of the SM Hierarchy problem. 
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Table 1: Tuning measures for models with only one parameter and one observable 



V Toy Models 

We now compare some of the tuning measures for various toy models and discuss the 
implications. In each of these examples we will assume a uniform probability distribution 
for the parameters. 

Table [T] compares the analytical results of various tuning measures for the simple mod- 
els with only one parameter and one observable. With only one parameter it is trivially 
the case that = A^g, so it is not included. 

In the first row of Table [T] are the results for a toy version of the Standard Model 
Hierarchy Problem, where we know the tuning is enormous. Here there is only one 
observable, the physical Higgs mass, mn- At one loop we write, 

m]j = ml- CA^ (18) 

and treat only the bare mass squared, rn^ as a parameter. A, the Ultra- Violet cutoff, is 
taken to be the Planck Mass or some other fixed scale, while C is a positive constant. 

Our measure was obtained by simply varying the tree-level mass parameter over the 
arbitrary range [ant^, hrn^ and applying the same dimensionless variations to the observ- 
able. This gives F = [b — a)ml and G = [b — a)m'jj, leading to the result for A shown. 
Notice that the arbitrary range [a, b] has fallen out of the result and it matches that 
obtained using the traditional measure, as shown earlier for all linear functions. 

We also determine A, and A^^. In both cases this introduces a dependence on the 
allowed range of ml in the theory, so we specify ^omm — — '^omax; present results 
where ml^-^ > CA^, though similar results can be obtained for other scenarios. These 
bounds give the total allowed range of the parameter in this model and should not be 
confused with the range of dimensionless variations which appears in the definition of F. If 
we take the range of variation to be large, ""^omaa; ~ ""^omm ^ CA"^, then A ^ ^ = Abg- 
Alternatively, if we choose a very narrow range of variation about CA^ + fi'j^, where 
fiH ~ 100 GeV, then A is very small. 

This is intuitively reasonable. Imagine some compelling theoretical reason for the 
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bare mass to be constrained close to the cutoff, e.g. CA^ < itlq < CA^ + (150 GeV)^. In 
light of this, the case for new physics at low energies would be dramatically weakened. 
Indeed it is precisely because there is no such compelling reason that we worry about the 
hierarchy problem and look to BSM physics such as supersymmetry to explain how we 
can have mn <^ ^pianck- 

Now let us compare this with the result for A^c- ^Hmm and m/^max are the extremum 
values of uih, dictated by the extremum values of mo. Notice that as mo^ax ~^ '^omm 
we have niHmin and mHmax — ^ f^H, so Aac 1- However a fundamental dif- 

ference between our measure and Aa_c is that the latter will give a large tuning for any 
mjj- -C K^^omaa; + "^Omm) ~ CA'^. If the Upper bound is chosen such that, ml^^^ > mg, 
then even a Higgs mass of 0{mQ) will appear fine tuned. This measure is not sensitive 
to the unnatural cancellation which causes our concern. Instead it is sensitive to the fact 
that large values of mfj take up a much larger volume of parameter space than small 
values of m^. This would be true even if the Higgs mass was described by mfj = rriQ, 
with no unnatural cancellation. 



Also shown in Tabled] are the results for the simple functions f{x) = x" and g{x) = 
e^^ . Earlier we showed there was no relative sensitivity in f{x). While Abg and A 
can be large despite the absence of relative sensitivity, our measure. A, is exactly unity 
for all X. Anderson and Castano's measure does remove the global sensitivity, but their 
tuning criterion prefers the observable to be as large as possible. For g{x), while there is 
relative sensitivity between different values of x, the constant factor k makes A^g* > 10 
for all |x| > 10/fc. For A the situation is similar, with A > 10 for all |x| > 10/i^, where 
K = k{b — a)/ln^. In A the effect of K is removed and though tuning still increases 
with X, this is now contextualised by comparing it to A. It is interesting that our mea- 
sure considers f{x) to have consistently no tuning (A = 1), whereas it is for g{x) that 
Aac = 1 for all x. 



The original illustration of global sensitivity presented by Anderson and Castano in 
Ref. f32] was for the proton mass. The proton can be much lighter than the Planck Mass 
without fine tuning because the renormalisation group equations (RGE) lead to only a 
logarithmic dependence on high scale quantities. However, by using the one loop RGE 
for the QCD coupling, 0:3, and equating the proton mass to the QCD scal^ 



Mproton ~ AqcD = C CXp 



(19) 



where 5^3 is the strong gauge coupling evaluated at the Planck scale, Mpianck, and C is a 
positive constant. As they demonstrated, this gives ABcigs) > 100. 

The analytical results for tuning in the mass of the proton, using Eq. (fT9|) . are shown 
in the final row of Table [TJ Notice that while the unnormalised tunings are both > 100, 
Aac and A are small. The latter has been determined only approximately in the limit 
where g^lnb <ti k and g^lna k for all gmin < (73 < g-max, where k = Sn'^/b^. 



In these one parameter examples the need for a normalised tuning measure is apparent. 
However A^c diverges significantly from our new measure, which in many of these simple 

iPor details see ^ 
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one dimensional models is equivalent to normalising the traditional measure with it's 
mean value. 

It is also interesting that even after accounting for global sensitivity some of these 
one dimensional functions may still show some small degree of tuning. This opens up the 
possibility that changing the parameterisation of the effective low energy theory might 
exacerbate or alleviate the tuning problem. Finding choices of parametrisation which 
reduce tuning could allow us to select high scale theories which are preferential in terms 
of naturalness. This point has not appeared in the literature and merits investigation. 
However we do not address this here but leave it for a future study. 

Now we consider models with more than one parameter. In these cases diverges 
from Abg and we must compare each of these with A. 

First we return to the SM hierarchy problem, but this time treat mn as a function 
of two parameters, m-Q and A^. In the one dimensional example the tension between the 
weakness of gravitation (the large Planck Mass) and a light Higgs mass was examined 
indirectly by choosing the Planck mass to be a fixed constant in theory. We now take a 
more direct route with two observables and Mp^^^^^. ("observed" to be large due to 
the weakness of gravitation), predicted from the parameters with, 

MLnc^ = A^ mj, = ml- CA'. (20) 

We are still predicting m'jj from Eq. (ITSll and have not split up any of the terms to 
introduce new cancellations, so we expect to simply reproduce the same result for A as 
we obtained in the one parameter toy SM model. However, the method applied provides 
a simple illustration of how our measure works with more than one parameter. We have 
a two dimensional parameter space, so allowing the parameters to vary about some point 
P'(mQ, A^) over the dimensionless interval [a, b] defines an area, F, in this space. Clearly 
the bounds from dimensionless variations in Mp^^^^^ are the same as those from A^, while 
the bounds from dimensionless variations in mjj introduce two new lines in the parameter 
space. 

This is shown in Fig. Illll for two different points. In the first point, the values of the 
parameters are of the same order as the observable, m^, because we have chosen a small 
value of Mpianck- So G is not much smaller than F. For the other point Mp^^^^^^ ^ 171^, 
resulting in an F much larger than G and fine tuning. Of course neither of these points are 
representative of the weakness of gravitation we observe. A point with Mpianck = 10^^ GeV 
and ruH = 120 GeV, would have F ^ G to such an extent that a graphical illustration is 
not possible. 

In general the areas are, F = (h — aYm^K^ and G = (h — aYh^m?^ so, 

A = 1 + ^ = /\bg. (21) 



H 



In this simple case we find the same result as the traditional measure. Combining Abg(A) 
and /XBGij^V) as if they are uncorrelated, gives. 



m: 



A^ = ^ (22) 
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Figure III: The two dimensional volumes (areas) F (dark grey) and G (light grey) for 
two different points in the two dimensional parameter space. 



With CA^ and itlq both ^ m\j, i.e. fine tuned scenarios, this gives us /\e ~ 72 A. While 
our measure does not deviate from /\bg in fhis simple example, models with additional 
parameters allow the observable to be obtained from cancellation of more than two terms, 
complicating the fine tuning picture. 



We now look at a model with four observables, M^, , M|, M|, and three parameters, 
Pl^phph described by, 

= cip\ - C2pI + CspI (23) 

M^=pI MI=pI MI=pI (24) 

For a point {m\, m\) , in the three dimensional parameter space, the traditional mea- 
sure gives AsaiPi) = Cimf/M'^ (no sum over i is implied), so. 



Abg = niax 



M2 



and Ae 



E 



2 4 



M2 



(25) 



To apply our tuning measure in the three dimensional case we must determine volumes 
F and G. For a point, (m^, mg, Tn"^), with = Mq = Cim^ — 02^2 + c^fn^ we have. 



d^F 



dpidpldpl 
d^G 



\[e{pl - arrrl)e{hml - p 



i=l 



d^F 



(26) 
(27) 



dpidpldpi dpfdpldpl 

where the latter uses Mf = pf and 6{x) is the usual Heaviside step function. Integrating 
Eq. (l26l) over all three pi gives the volume. 



(28) 



and similarly Eq. f l271) gives, 

G = {h — a)^ < 9{c'ira\ — C2ra\)9{c2ra\ — Cim\) 



C3 



-m\mlM'^ 



3C2C3 
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+ e{cim\ - C2ml)e{c2ml - cml) 
+ e{cW^ - C2ml)e{cmi - c^ml) 
+ e{c2ml - c,ml)9{c2ml - cW^) 



3C2Ci 



-rrin 



2 2 2 

mim2mr^ 



3CiC3 



-■rrin 



^ 2 2 (1 /r2 ^ 

— mim^M 

C2 3C1C2C3 



/r6 



(29) 



We find that the analytical expressions for tuning in this model depend on the mass 
hierarchy of mi, 1712 and m^. 



For Ciml > 02^^ > c^m'^ we find, 



A = S 



For C3m| > 02^2 > Ciml we find: 
A 



Abg if csmg -C C2m2. 



(30) 



A 



BG 



if 



(31) 



For CsmI > Cimf > C2Tn\ and Cim\ > c^m\ > C2Tn\. 



A 



2 2 
mim3 



2 2 



■( 2 2^-^2 4 

if Cim^c^rrij^ ^ C2''^2- 



(32) 



For 02^2 > Cim^ > 037723 and C2m2 > 037773 > Cim\: 



A 



C2'rn\m\m\ 



G m'im%M^ - 



Abg if -C C17771C37773. 



(33) 



Notice that these results do not match A^;, but in three dimensions at least A^g is a 
much better approximation, as is shown in Fig. IIVI 

However, as we have seen, in moving from two parameters to three parameters these 
discrepancies appeared, increasing the number of parameters further will increase the 
divergences between the measures. 



VI Fine Tuning in the MSSM 

The analytical methods described above become increasingly complicated to apply as the 
number of parameters and observables are increased. For such situations we have also 
developed a numerical procedure which can be applied to produce approximate results 
for tuning. Since the MSSM contains many parameters and many observables we chose 
to apply our numerical approach here. 

We take random dimensionless fluctuations about an MSSM point at the GUT scale, 
P' = {Pk}, to give new points {Pi}. These are passed to a modified version of Softsusy 
2.0.5[48j. Each random point Pi is run down from the GUT scale until Electroweak 
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Figure IV: Comparison of (unnormalised) tuning measures in the three parameter model 
with m| varying from to 100 GeV^ and Mq = 1 GeV^ and = 99 GeV^ kept constant. 
ml then varies according to Eq. ([23]) to accommodate the changes in mg. Left: between 
Abg and our new measure. Right between A^; and our new measure. 



Symmetry is broken. An iterative procedure is used to predict M| and then all the 
sparticle and Higgs masses are determined. For a theoretical discussion, see Ref . [19] • 

As before F is the volume formed by dimensionless variations in the parameters. Go, 
is the sub-volume of F additionally restricted by dimensionless variations in the single 
observable Oi, a < 'o\{l/}] — ^- -^^ usual G is the volume restricted by a < Q^^j^f |^ < b, 
for each observable, Oj, where {Oj} is the set of masses predicted in Softsusy. For every 
Oi a count, iVo., is kept of how often the point lies in the volume Gpi as well as an overall 
count. No, kept of how many points are in G. Tuning is then measured according to, 

Ao. ^ (34) 



for individual observables and 



for the overall tuning at that point. 



A « - (35) 



Before describing the results two comments on this approach should be made. Firstly 
when using Softsusy to predict the masses for the random points, sometimes problems 
are encountered. We may have a tachyon, the Higgs potential unbounded from below, or 
non-perturbativity. Such points don't belong in volume G as they will give dramatically 
different physics. However it is unclear which volumes. Go,? the point lies in. Such points 
never register as hits in any of the Go, and this may artificially inflate the individual 
tunings, including 1^m\- Keeping the range small reduces the number of problem points. 
Therefore we chose a = 0.9 and b = 1.1 for our dimensionless variations. 

Secondly, since we are measuring tuning for individual points numerically and cover 
only a small sample of points, it is not possible to obtain mean values of A and the Ao, 
as we haven't sampled the entire space. When simply comparing how the tuning varies 
about the parameter space the normalisation factor is not needed. However to compare 
the tuning between different observables as well as to compare with different models some 
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form of normalisation is essential. 



We considered points on the Constrained Minimal Supersymmetric Standard Model 
(CMSSM) benchmark slopfl SPS la [SU]. This slope is defined by, 

mo = —Aq = 0.4mi, sign(yu) = +, tan (3 = 10, (36) 

where mo is the common scalar mass, mi/2 the common gaugino mass (both at the GUT 
scale) and sign(/i) is the undetermined sign of fi, the magnitude being determined from 
a loop corrected, inverted form of Eq. ([T]) with M| set to it's observed value. Aq is 
the common multiplicative factor which relates the supersymmetry breaking matrices of 
trilinear mass couplings to their corresponding Yukawa matrix, e.g. = Aqi/u. 

The parameters we vary simultaneously are the selEl {mo, mi/2, I^gut, Aq, yt, y^, y^-}, 
where is the soft bilinear Higgs mixing parameter and yt,yb, Ut are the Yukawa cou- 
plings of the top, bottom and tau respectively. The gauge couplings are not included as 
parameters. Doing so would introduce excessive global sensitivity, increasing the statistics 
needed to keep the errors under control. 

First we applied our tuning measure to the observable M| for 13 points on the SPS la 
slope. Moving along this slope in mi/2 is an increase in the overall supersymmetry 
breaking scale, since the magnitude of every soft breaking term is increasing. We have 
plotted the results of this investigation in Figure |Vl 



500 
400 
300 

z 

200 
100 



200 400 600 800 1000 

m,,JGeV) 




Figure V: for the SPS la slope. Error bars denote a one standard deviation statistical 
error arising from the numerical procedure. 

As expected there is a clear increase in tuning as the supersymmetry breaking scale is 
raised. The statistical error also increases with the tuning, making the numerical approach 

^Such benchmark slopes and points, known as Snowmass Points and slopes (SPS)[SD] are chosen by 
consensus as representing qualitatively different MSSM scenarios and are very useful for comparison with 
other work. 

•^Note that since points on the SPS la slope have \fj,\ set by M^, our tuning measure is not sensitive 
to the /i-problcm. However for our random variation about the SPS la points we do treat (J-gut as a 
parameter because we are predicting AI^ from the parameters, not fixing it to it's observed value. 
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most difficult to apply when the tuning is large. However precise determinations of tuning 
are only relevant for moderate and low tunings. With tunings greater than 500, precise 
values are not required. 

Due to the difficulty in this approach for measuring large tunings we looked in more 
detail at points expected to have moderate tuning. We chose a grid of points with, 

Ao = -100GeV, tan/3 = 10, sign(/i) = +, 
250 GeV < m 1 < 500 GeV, 100 GeV < niQ < 200 GeV. (37) 

Shown in Fig. IVII (top) is a plot of Ajvj2 over this grid of points. While the errors are 
still significant (< 10%) there is a clear trend of tuning increasing with mi/2- Also shown 
(bottom left) is Aj^2 averaged over the five different values of niQ. This substantially 
reduces the errors giving a much more stable picture of tuning increasing linearly with 
mi/2- Similarly A^.^i , averaged over the eleven different values of mi/2, is shown (bottom 
right) as a function of mg. appears insensitive to variations in ttiq. These trends 

can be understood by looking at the one loop renormalisation group improved version of 
Eq. ([T]), written in terms of the the parameters (with tan (3 = 10), 

M| ^ 2(-|/i|2 + 0.076m^ + 1.97mi + O.IOA^ + 0.38Aomi ), (38) 

2 2 

where is the value at Mz and and differs from the parameter at the GUT scale, 
fiGUT- The large coefficient in front of mi/2 explains why explains why variations in this 
parameter have a much greater impact on Ajv,/2 than variations in mo whose coefficient 
is much smaller. 

A, which includes all of the masses predicted by Softsusy as well as M|, is shown 
in Fig. IVIII Although the errors are much larger here, a similar pattern to that for 
M| can be seen. Since these are unnormalised tunings, the numerical values of the two 
measures cannot be compared and one should not assume that A > A^j2 implies that 
the tuning is worse than when only M| was considered. In fact the lack of evidence for 
distinct patterns of variation in tuning from the Figs. |Vl] and IVIII is consistent with the 
conjecture that the large cancellation between parameters in M| is the dominant source 
of the tuning for these points. 

Fig. IVIIII shows that A^2 and A^2 have similar patterns of variation to Aj^.f2 and 

A over mi/2, though the gradients are noticeably shallower. While we know m\ and m?_^ 
contribute to the Little Hierarchy Problem by giving a large contribution to M|, thereby 
requiring a cancellation to keep Mz light, this shows there is also some tension in their 
own masses which restricts the parameter space. It is not clear from our results whether 
or not dimensionless variations are restricting different regions of parameter space to 
those in M| or if G^2 and G^2 are merely sub- volumes of (^^^2 , with no influence on 

A. This topic deserves further study. 

However our results do show some evidence that the Little Hierarchy Problem is 
not the only source of tuning. Displayed in Fig. [IX] is Ajy^2. Notice that A^j2 is very 
small, so the errors are significantly reduced and we can resolve very small variations in 
Ajv^^. As with the other observables tuning increases with mi/2, but it is a distinctly 
non- linear variation. More surprising is that tuning decreases with mo. This pattern of 
variation, distinct from that shown for A]^2 , shows a different source of tension. It can 
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Figure VI: Tuning variation in M|. Top: Aj^fZ for all points on our grid. Bottom left: 
Ajy^l plotted against mi/2- To reduce statistical errors, at each value of mi/2, we have 
taken the mean value A^^i over the five different ttlq values. Bottom right: Aj^/2 plotted 
against uiq. To reduce statistical errors, at each value of mo, we have taken the mean 
value Aj^,/2 over the eleven different mi/2 values. 

be understood by examining the one loop RGE solution for M^, 

Ml ^ 2f{\^iGUT\\ [gi], {Vi]) + 0.81m^ - l.bbml - 0.022/1^ - 0.41Aomi , (39) 

2 2 

where / is a function of supersymmetry preserving parameters only, arising from the 
evolution of Notice that there is some opportunity for a cancellation here to make 
Ma lighter than expected. However the cancellation in the points we have looked at 
is very small, leading to small values for /^m\- As increases the already dominant 
positive part of the equation increases and M^ increases. As this happens the cancellation 
becomes less significant to Ma further reducing A^2 as shown in Fig. IIX( bottom right). 
Increasing mi/2 increases the size of the cancellation. If all other parameters on the right 
hand side of Eq. (l39l) were fixed then we would expect to see A^,/2 increase linearljfl with 
mi/2- However each point on our grid has the value of Mz = 91.188 GeV fixed, and the 

''The effect of Aomi/2 can be neglected since mi/2 > A^. 
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Figure VII: Variation in A plotted as in Fig. IVII for 
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Figure VIII: Variation of unnormalised tunings in the mass of the heaviest stop (A 
shown left) and the mass of the light Higgs (A^2, shown right) over mi/2 
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term f{\^\^,{gi},{yi}) ^ changes according to an inverted Eq. ( 138|) . This means 
M\ is also increasing with mi/2 and the balancing act between these two different effects 
leads to the nonlinear pattern shown. 
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Figure IX: Variation in Ajyj2 plotted as in Fig. |Vl] for A 



Although we can't determine the normalisation using this approach it is nonetheless 
interesting to compare the unnormalised tunings for the points in our study with those 
obtained for points with more "natural" looking spectra. We present two points for this 
purpose. NPl and NP2 are defined by, 



NPl : mi 
NP2 : ml 

2 



-Mz 



mo = Mz, ao- 
50GeV, mo = 100GeV, ao- 



-Mz, 



sign(//) = +, tan/3 = 3, 
-50GeV, sign(/i) = +, tan/5 = 10. 



(40) 



The spectra of these points are displayed in Fig. |X] and Fig. IXIl and the unnormalised 
tunings are displayed in Table O Note that these are not intended to be "realistic" 
scenarios. Indeed both NPl and NP2 are ruled out by experiment but are simply intended 
to provide "natural" scenarios for comparison. 

While NPl has low values of A»^2 A™2 A™2 and A»^2 it has a relatively large 

Z h i2 ^ 

tuning in the mass of the lightest neutralino (A^2 ). These combine to give a A which 

^1 

is similar in size to the values found for our grid of points. In NP2 all of the tunings 
are relatively small, but the combined tuning is still larger than may naively have been 
anticipated. This is because many of these small tunings for individual observables are 
not correlated and are restricting different regions of parameter space. Table [3] shows the 
approximate relative magnitude of the tunings in our grid points with respect to these 
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Figure X: Point NPl with a "natural" spectrum 
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Table 2: Unnormalised tunings for the two points, NPl and NP2, with natural looking 
spectra. 



seemingly natural points. 

In attempts to find a CMSSM scenario with a mass spectrum which is manifestly 
natural we found many scenarios where tuning appeared in the mass of the lightest 
neutralino. NPl is a (moderate) example of this. This is because in some parameter 
choices, the lightest neutralino becomes very light due to large cancellations between 
the parameters. Other observables may also contain large cancellations between the 
parameters in certain regions of parameter space. While we have not studied this enough 
to make definitive claims, this may suggest that mass hierarchies appear in a greater 
proportion of the parameter space than conventional CMSSM wisdom dictates. This 
would reduce the true tuning in the CMSSM as scenarios with hierarchies would be less 
atypical than previously thought. A reduction in tuning from this effect can only be 
measured by using our normalised new measure, A. 

Unfortunately the numerical approach we have applied to the MSSM in this paper 
cannot be used to address this issue. An average measure of A, over the whole parameter 
space, is needed in order to investigate this possibility. A thorough numerical survey 
of the parameter space would be too expensive, however an analytical study may be 
more promising. Findings in numerical studies like this may be used to identify which 
observables and parameters are important for fine tuning and therefore reduce the set 
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Figure XI: Point NP2 with a "natural" spectrum 





A 


Ami 


A„^-^ 


A 2 

n. 


Ami 


A 2 


Relative to NPl 


0.5. .1.5 


3.. 10 


1..2 


7.. 25 


1 


0.2 


Relative to NP2 


5.. 15 


10. .50 


4.. 7 


6. .23 


1 


2 



Table 3: Approximate relative tunings for the points in our study, with respect to those 
for NPl and NP2. 



{Oi} and {pi} to a manageable size. We will not carry out this programme here, but 
leave it for a future study. 

It is not just the possibility of finding a larger than expected global sensitivity which 
motivates this study. It may be that most of the CMSSM parameter space is hierarchy 
free and this is not a significant effect. However identifying a region of parameter space 
where mass hierarchies are common also opens up new possibilities. Past studies (see 
e.g. [5T|52] ) have looked for a theoretical basis for relations between parameters which 
enforce a hierarchy between Mz and Msusy- However no search has been made for the- 
oretical relations which simply restrict the parameter space to regions where hierarchies, 
in general, are common. Such studies may also have the possibility of solving the Little 
Hierarchy Problem. 

Here we have two complimentary approaches. An analytical approach which can de- 
termine tuning precisely, but is complicated and unwieldy when applied to a great number 
of parameters and observables and a numerical approach which can be applied to such 
situations but is not able to give an unambiguous measure of tuning as global sensitivity 
cannot be accounted for. Progress can be made by combing our two approaches. Since 
solving for the tuning analytically with all parameters and observables included would 
be difficult, one should first apply the numerical method. This might identify which ob- 
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servables are in tension and responsible for the restriction of parameter space and also 
along which axis in parameter space this restriction takes place. If these are a sufficiently 
small set (maybe no more than 5 parameters and 5 observables) then the analytical mea- 
sure can be applied to this limited set to obtain a reasonably accurate and unambiguous 
measure of tuning for that model. 



VII Conclusions 

Fine tuning ^ 10^^ within the Standard Model has motivated many of the BSM the- 
ories which are popular within particle physics. In particular it motivates low energy 
super symmetry. However constraints from LEP and other searches have placed stringent 
bounds on new physics which mean that many of the proposed solutions to the SM fine 
tuning problem also require tuning to some degree. In order to compare the viability of 
such models and judge whether or not they are satisfactory a reliable measure of tuning 
is required. 

Current measures of tuning have several limitations. They neglect the many param- 
eter nature of fine tuning, ignore additional tunings in other observables, consider local 
stabihty only and assume Csusy is parametrised in the same way as Cqut- In the hter- 
ature there have been different approaches to combine tunings for individual parameters 
and observables. With no guiding principle to select one particular approach, which 
models are preferred in terms of naturalness can depend on which tuning measure is 
used. 

In this paper we have presented a new measure of tuning based upon our intuitive 
notion of the restriction of parameter space. This measure can also be obtained by gener- 
alising the traditional measure of tuning to include many parameters, many observables 
and finite variations in the parameters followed by removing global sensitivity by factoring 
out the mean value of the unnormalised sensitivity. 

From the application of this new measure to various toy models, we have shown 
that none of the other measures satisfactorily combine individual tunings per parameter. 
Interestingly though, in the absence of global sensitivity, it is the traditional measure of 
Barbieri and Guidice which comes closest to our result with deviations for these simple 
examples being very small. 

A numerical approach for some CMSSM scenarios demonstrated how the tuning in 
complicated models with many parameters and many observables may be examined and 
also highlighted some of the complications and issues encountered in doing so. 

Our new measure is needed in future studies to examine tuning in the Z boson mass 
and cosmological relic density simultaneously; to judge the true tuning in the NMSSM 
in light of to examine parametrisation choices which alleviate the tuning in different 
models and to study the global sensitivity of the complete tuning measure to see if this 
may cause a significant reduction in the tuning problem. 
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